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ABSTRACT 



We examine the excitation of oscillations in the magnetic network of the Sun through 
the footpoint motion of photospheric magnetic flux tubes located in intergranular lanes. 
The motion is derived from a time series of high-resolution G band and continuum 
filtergrams using an object-tracking technique. We model the response of the flux 
tube to the footpoint motion in terms of the Klein-Gordon equation, which is solved 
analytically as an initial value problem for transverse (kink) waves. We compute the 
wave energy flux in upward propagating transverse waves. In general we find that 
the injection of energy into the chromosphere occurs in short-duration pulses, which 
would lead to a time variability in chromospheric emission that is incompatible with 
observations. Therefore, we consider the effects of turbulent convective flows on flux 
tubes in intergranular lanes. The turbulent flows are simulated by adding high-frequency 
motions (periods 5-50 s) with an amplitude of 1 km s _1 . The latter are simulated by 
adding random velocity fluctuations to the observationally determined velocities. In 
this case we find that the energy flux is much less intermittent and can in principle 
carry adequate energy for chromospheric heating. 

Subject headings: MHD — Sun: magnetic fields Sun: oscillations 



1. INTRODUCTION 

It is now well established that the solar photosphere is threaded with strong (kilogauss) mag- 
netic fields in the form of vertical flux tubes located in the network, where they are observed as 
bright points. Observations have revealed that these network bright points (NBPs) are in a highly 
dynamical state due to the buffeting action of granules (e.g. Muller 1983; Muller et al. 1994; Berger 
& Title 1996; Berger et al. 1998; van Ballegooijen et al. 1998). 

The jostling of magnetic elements by granules can excite oscillations in flux tubes, which can 
in principle be an important mechanism for chromospheric and coronal heating, provided that 
the motions are rapid enough (Choudhuri, Auffret & Priest 1993; Steiner, Knolker &; Schiissler 
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1994; Hasan & Kalkofen 1999, hereafter Paper I). As transverse waves, excited in the photosphere, 
travel upwards their velocity amplitude increases. In chromospheric layers, the velocity amplitude 
becomes comparable to the tube speed for kink waves, leading to an efficient coupling with longi- 
tudinal waves (Kalkofen 1997). The latter can dissipate by forming shocks (e.g. Zhugzhda, Bromm 
& Ulmschneider 1995). 

In Paper I it was argued that network oscillations can be efficiently excited through the buf- 
feting action on flux tubes by external granules with intermittent motions (of about 2 km s _1 ) 
occurring on a time scale of less than about half the cutoff period of kink waves. Typically, a pe- 
riod of 1 min. is consistent with the general result. We assumed an analytic form for the external 
motions. In the present work, we use observational data obtained at the Swedish Solar Observatory 
at La Palma during 1995 (Berger et al. 1998) to infer the horizontal footpoint motions of magnetic 
elements in the photosphere. We model the response of flux tubes to these motions and calculate 
the energy flux in upward-propagating transverse (kink) waves. We also simulate the effect of hy- 
pothetical motions on time scales of a few seconds, which could significantly enhance the energy 
flux in kink waves. Finally, we examine some implications of our results. 

2. MEASUREMENT OF BRIGHT POINT MOTION 

The data were collected on October 5, 1995 between 10:57 and 12:08 UT, simultaneously in 
the G band (4305 A) and the nearby continuum (4686 A), using the phase-diversity method (see 
Lofdahl et al. 1998; Berger et al. 1998). After filtering out p— mode oscillations, a time series 
of 180 images was obtained, with high spatial resolution and with an interval between frames of 
23.5 s. Following Berger et al., we subtract the images in the G band and 4686 A, which clearly 
highlights the network bright points. 

The motion of the bright points was followed using a tracking technique with "corks" as tracers 
of bright points (van Ballegooijen et al. 1998, hereafter Paper II). The corks, each with a radius 
of 60 km (the pixel size of the current data set), were advected by an artificial flow field that 
is proportional to the gradient of the intensity in the magnetic image, and their positions were 
followed in time (for details see Paper II). In Figure 1 we show the cork positions superposed on a 
typical frame of a magnetic difference image, where x and y denote the horizontal coordinates on 
the Sun. Figures 2a and 2b show the x and y components of the velocity as a function of time for 
two representative corks, whose initial locations are indicated by arrows in Figure 1. 

3. KINK WAVE EXCITATION DUE TO FOOTPOINT MOTION 

Let us assume that the motion of the G band bright points can be taken as a proxy for the 
footpoint motion of flux tubes at the base of the photosphere (z = 0). Furthermore, for simplicity, 
we assume that flux tubes are isothermal and "thin" compared to the pressure scale height H. It is 
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convenient to work in terms of the "reduced" transverse displacement Q±(z,t) = ^±e~ z ^ 4H , where 
is the Lagrangian displacement. Once the footpoint motion is specified, the velocity (Q±) at 
any height z and time t can be determined by the following expression (Paper I): 




where J\ is a Bessel function, ( K = y/(t — t ) 2 — (z/c K ) 2 , c K and oj k are, respectively, the propagation 
speed and cutoff frequency for kink waves, given by: 




7 1 + 2/3 




_9 1_. 

8H 1 + 2/3' 



c s is the sound speed, k K = uj k /c k , 7 = 5/3, /3 = Sirp/B 2 , p(z) is the gas pressure inside the tube 
and B(z) is the magnitude of the vertical component of the magnetic field on the tube axis. The 
energy flux (F K ) in a single flux tube is given by (Paper I): 



where p e $ denotes the external gas pressure at z = 0. 

Figures 3a and 3b show the vertical energy flux in transverse waves versus time at a height 
z = 750 km for two observed magnetic elements located initially at the coordinates shown above 
the figure. The left and right panels correspond to footpoint motions with the x and y components 
of the velocity, respectively. 

We find that the injection of energy into the chromosphere takes place in brief and intermittent 
bursts, lasting typically 30 s, separated by longer periods (longer than the time scale for radiative 
losses in the chromosphere) with lower energy flux. The peak energy flux into the chromosphere 
is as high as ~ 10 9 erg cm -2 s~ 1 in a single flux tube, although the time-averaged flux is ~ fO 8 
erg cm" 2 s _1 . This scenario for heating the magnetic network would produce strongly intermittent 
chromospheric emission consisting of brief intense flashes superimposed on generally a very low 
background. However, in dense network regions the observed chromospheric emission is stably 
present and exhibits relatively low amplitude, long-period variations (e.g. Dame, Gouttebroze and 
Malherbe 1984; von Uexkiill et al. 1989; Lites, Rutten and Kalkofen 1993). Therefore, the above 
scenario cannot readily explain the observed persistence of emission from dense network regions. 
A possible remedy is to consider the effect of high frequency motions. 
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4. SIMULATING HIGH-FREQUENCY MOTIONS 

The observations of G band bright points were taken with a time cadence of 23.5 s. However, 
the possibility of random footpoint motions occurring on a shorter time scale cannot be ruled 
out. The solar convection has very high Reynolds number and is expected to be highly turbulent, 
involving a wide range of length and time scales. Indeed, spectroscopic observations of the solar 
granulation have shown enhanced broadening of spectral lines in intergranular lanes, indicating 
enhanced turbulence within these region (e.g. Nesis et al. 1992, 1994). Magnetic flux tubes interact 
with these turbulent flows within the lanes. We speculate that these interactions produce transverse 
displacements of magnetic field lines on length scales much less than the flux tube diameter (~ 100 
km), i.e., the magnetic field inside a flux tube has fine structure on transverse scales £ <C 100 km. 
These perturbations likely generate Alfven waves which propagate upward along the field lines and 
dissipate their energy higher up in chromosphere. 

The velocity at small spatial scale can be estimated as follows. Assuming the turbulence has 
a Kolmogoroff spectrum, the velocity of eddies on length scale I is given by v(£) ~ vq{1/ Iq) 1 ^ , 
where £q is the outer scale of the turbulence, and vq is the velocity of the largest eddies. We assume 
that £q is given by the width of an intergranular lane as determined from 3D simulations of the 
solar granulation (e.g., Stein and Nordlund 1998), and that vq is the convective velocity predicted 
by such simulations: £q ~ 100 km and vq ~ 2.5 km s -1 . This turbulence model predicts that 
v ~ 1 km s _1 on a length scale £ ~ 10 km, i.e., there is significant power in motions at this scale. 
The time scale of such motions is £/v ~ 10 s, which is significantly shorter than the time scale for 
granular buffeting. 

In the following we simulate the effect of such high-frequency motions on the response of the 
flux tube (similar to the study of longitudinal wave excitation in flux tubes by Ulmschneider and 
Musielak 1998). For simplicity we ignore the small transverse scale of the motions and we describe 
the "turbulence" in terms of footpoint displacements of the flux tube. We do this by superposing 
random motions, with a rms amplitude of 1 km s _1 and a zero mean value, on the observationally 
determined velocity. The random motions vary on a time scale of 2.35 s. Such high-frequency 
motions (if present) would not have been detected in the present observations. Hence, it is not 
possible a priori to rule them out. Results from such simulations are shown in Figures 4a and 4b 
for the x component of the velocity as a function of time for two representative corks whose initial 
locations are indicated above each panel, and in Figures 5a and 5b for the corresponding vertical 
energy wave flux. We note that the energy flux in upward-propagating transverse waves is much 
larger and less intermittent than in the case without high-frequency footpoint motions. 

The energy flux in the form of high-frequency transverse waves appears adequate to contin- 
uously heat the chromosphere. Taking a filling factor of 10% at a height z = 750 km above the 
photosphere, the predicted flux is approximately 2 x 10 7 erg cm~ 2 s _1 , which is sufficient to balance 
the observed radiative loss of the chromospheric network (~ 10 7 erg cm -2 s _1 , see Model F' in 
Avrett 1985) 
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5. DISCUSSION 

The energy flux in upward propagating transverse waves has been shown for a few represen- 
tative magnetic elements. We should point out that we carried out similar calculations for several 
corks, before selecting two as being typical of the large sample that we examined. 

We find that the injection of energy into the chromosphere occurs in brief and intermittent 
bursts (lasting typically 30-60 s), separated by longer periods (longer than the typical time scale 
for radiative loss) with low energy flux. The peak energy flux into the chromosphere is typically 
10 9 erg cm~ 2 s _1 in a single flux tube. Observationally, such a scenario for heating of the magnetic 
network on a 1 arc sec scale would lead to high variability in Ca II K emission, contrary to 
observations. However, the observations of G band bright points used here do not have the time 
resolution and sensitivity necessary to detect motions on shorter time scales. We simulate such 
high-frequency motions by adding random velocity fluctuations with a zero mean value and an 
rms amplitude of 1 km s^ 1 . We conclude that for transverse waves to provide a viable mechanism 
for sustained chromospheric heating, the main contribution to the heating must come from high- 
frequency motions (periods 5-50 s). 

The proposed high-frequency motions have very small spatial scale (10 - 100 km) and the 
transverse displacements involved are also very small (for example, a motion with period 50 s and 
velocity amplitude of 1 km s _1 has a displacement amplitude of only 8 km). Such length scales 
are well below the angular resolution of present-day solar telescopes. Therefore, it is unlikely that 
such high-frequency motions will ever be detected in studies of the proper motions of photospheric 
flux tubes. Spectroscopic methods, such as measurements of spectral line broadening by unresolved 
turbulent flows, offer greater promise but will still require very high angular resolution. 

The observational data used in this study were provided by R.A. Shine of Lockheed-Martin 
and were obtained at the Swedish Vacuum Solar Telescope. We acknowledge support from the 
National Science Foundation and NASA. SSH is thankful to the Smithsonian Institution for travel 
support under SMIN grant no. 9411FR00007. We are grateful to the referee for helpful comments. 



-6- 



REFERENCES 

Avrett, E. H. 1985, in Chromospheric Diagnostics and Modeling, ed. B.W. Lites, Nat. Solar Obs. 
(Sunspot, N.M.), p.67 

Berger, T. E., Lofdahl, M. G., Shine, R. A., & Title, A. M. 1998, ApJ 495, 973 

Berger, T. E., & Title, A. M. 1996, ApJ 463, 365 

Choudhuri, A. R., Auffret, H. & Priest, E. R. 1993, Sol. Phys. 143, 49 

Dame, L., Gouttebroze, P., and Malherbe, J.-M. 1984, A&A, 130, 331 

Hasan, S. S. & Kalkofen, W. 1999, ApJ 519, 899 (Paper I) 

Kalkofen, W. 1997, ApJ 486, L145 

Lites, B. W., Rutten, R. J., and Kalkofen, W. 1993, ApJ, 414, 345 

Lofdahl, M. G., Berger, T. E., Shine, R. A., & Title, A. M. 1998, ApJ 495, 965 

Muller, R. 1985, Sol. Phys. 85, 113 

Muller, R., Roudier, Th., Vigneau, J. & Auffret, H. 1994, A& A 283, 232 

Nesis, A., Hanslmeier, A., Hammer, R., Komm, R., Mattig, W., and Staiger, J. 1993, A&A, 279, 
599 

Nesis, A., Hammer, R., Hanslmeier, A., Schleicher, H., Sigwarth, M., and Staiger, J. 1996, A&A, 
310, 973 

Stein, R. F., and Nordhmd, A. 1998, ApJ, 499, 914 

Steiner, O., Grossmann-Doerth, U., Knolker, M., & Schiissler, M. 1998, ApJ 495, 468 
Ulmschneider, P. & Musielak, Z. E. 1998, A&A 338,311 

van Ballegooijen, A., Nisenson, P., Noyes, R. W., Lofdahl, M. G., Stein, R. F., Nordhmd, A., 
Krishnakumar, V. 1998, ApJ 509, 435 

von Uexkiill, M., Kneer, F., Malherbe, J.-M., and Mein, P. 1989, A&A, 208, 290 

Zhugzhda, Y. D., Bromm, V. & Ulmschneider, P. 1995, A&A 300, 302 



This preprint was prepared with the AAS IATjrjX macros v5.0. 




5 10 15 20 

X (Mm) 



Fig. 1. — Difference between G band and continuum images in the initial frame. The black dots 
depict the corks, used for tracking the bright points. The arrows denote the locations of two 
representative corks. 
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Fig. 2. — Variation with time of the x and y components of the horizontal velocity of G band bright 
points. The coordinates (in Mm) of the corks at the initial instant (shown in Figure 1) are given 
above each panel. 
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Fig. 3. — Variation with time of the vertical energy flux in transverse waves at z = 750 km 
due to footpoint motions corresponding to initial locations (a) (18.7 Mm, 11. 9 Mm) and (b) 
(14.5 Mm, 11.0 Mm). The left and right panels correspond to motions associated with the x and y 
corrmonents of the velncitv. resnectivelv. 
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Fig. 4. — Variation with time of the horizontal velocity of G band bright points with random noise 
added for footpoint motions associated with corks initially at (a) (18.7,11.9) Mm and (b) (14.5,11.0) 
Mm. 
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Fig. 5. — Variation with time of the vertical energy flux in transverse waves at z = 750 km due to 
footpoint motions with random noise associated with corks initially at (a) (18.7,11.9) Mm and (b) 
(14.5,11.0) Mm. 



